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An alternative parameterization of it-matrix theory is presented which is mathematicaUy equiv- 
alent to the standard approach, but possesses features which simplify the fitting of experimental 
data. In particular there are no level shifts and no boundary-condition constants which allows the 
positions and partial widths of an arbitrary number of levels to be easily fixed in an analysis. These 
alternative parameters can be converted to standard it-matrix parameters by a straightforward 
matrix diagonalization procedure. In addition it is possible to express the collision matrix directly 
in terms of the alternative parameters. 
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PACS numbers: 24.10.-i, 24.30.-v, 31.15.-p, 02.10.Ud 



I. INTRODUCTION 

The il-matrix theory of reactions has proven over the 
course of time to be very useful in nuclear and atomic 
physics, both for the fitting of experimental data and 
as a tool for theoretical calculations. In this paper we 
explore a mathematically-equivalent alternative formula- 
tion of i?-matrix theory which will be especially useful 
for the fitting of experimental nuclear physics data. 

In a recent paper an alternative parameterization R- 
matrix theory was described by Angulo and Descouve- 
mont . In their framework there are no level shifts and 
it is straightforward to incorporate known information 
about level energies and partial widths. They presented 
an approximate iterative relation between the alternative 
parameters and the standard i?-matrix parameters. In 
addition consideration was limited to the single-channel 
case with a boundary condition constant of zero. Some 
aspects of these alternative parameters have also been 
discussed in a paper by Barker [Q. In this paper we 
further develop the concept of an alternative i2-matrix 
parameterization. The description is generalized to allow 
non-zero boundary condition constants and an arbitrary 
number of channels. We present an exact method for 
converting the alternative parameters to the standard R- 
matrix parameters which only requires a matrix diagonal- 
ization. We also found a rather surprising result, that the 
collision matrix can be calculated directly from the alter- 
native parameters using alternative formulations of the 
level matrix or R matrix. We then discuss the solution of 
the nonlinear eigenvalue equation required to extract the 
alternative parameters from the standard parameteriza- 
tion, and demonstrate some of these ideas using a simple 
example. Finally we briefly discuss the application of the 
alternative parameterization to 7 rays and /3 decays. 



II. REVIEW OF STANDARD it-MATRIX 
THEORY 

We begin by reviewing some of the notation and results 
of standard i?-matrix theory as described by Lane and 
Thomas (LT) [pj. The R matrix is a function of the 



energy E and is defined by 

P _ 7Ac'7Ac 



(1) 



where E\ are the level energies, j\c are the reduced width 
amplitudes, A is the level label, and c is the channel label. 
We will assume that the numbers of levels and channels 
are finite and given by N\ and Nc, respectively. One 
must also specify the constants Be, which determine the 
boundary conditions satisfied by the underlying eigen- 
functions. 

In order to calculate physical observables one must 
employ various combinations of the Coulomb wavefunc- 
tions, evaluated at the channel radius Tc = Oc. The 
quantities Ic and Oc are defined by (LT, Eq. II. 2. 13). 
For closed channels the outgoing solution Oc is taken to 
be the exponentially-decaying Whittaker function (LT, 
Eq. II.2.17). In addition one defines fl^ = {Ic/Oc)'^/'^ 
and 



do, 



(2) 



where the shift factor Sc and penetration factor Pc are 
real quantities. The collision matrix U is an Nc x Nc 
matrix which determines the observable quantities; it is 
related to the R matrix via (LT, Eq. VII. 1): 

U = 2ip^l^O-^[l - R{L - B)]-'^Rp^l^O-^ + 10^ 

(3) 

where O, /, L, B, and p are purely diagonal with ele- 
ments Oc, Ic, Lc, Be, and kcac, respectively; 1 is the unit 
matrix, and kc is the wavenumber. 

It is convenient to form the level-space column vector 
'jc from the jxc, and to then form the rectangular matrix 
7 from the 7c such that the matrix 7 has N\ rows and 
Nc columns. In addition, the diagonal matrix e is defined 
by 

exf, = Ex5x^- (4) 

The R matrix defined by Eq. (Q) can now be written 
succinctly as 



R 



(5) 
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The collision matrix can also be expressed as 

U = 2ipl/2o-l^T^^^l/2(5-l _^ jQ 1 ^g^) 

where A is an N\ x A^^, matrix defined by its inverse: 

A-^ =e-El-j{L-B)j^. (7) 

The equivalence of these two forms for the collision ma- 
trix is discussed in (LT, Sec. IX. 1) and in the Appendix. 
In addition the elements of the collision matrix connect- 
ing open channels in Eq. ^ can also be expressed as 



Ur'r. — ^r.'^n 



(8) 



using the definitions of the Coulomb functions. 

An interesting feature of ii-matrix theory is that the 
collision matrix is invariant under changes in the Be, pro- 
vided that the E\ and 7ac are suitably adjusted. This 
result remains true even for the case of finite Nx [Q. 
The transformation is most easily described using ma- 
trix equations in level space. Let us consider the trans- 
formation Be ^ B'^, Ex E'^, and 7ac I'xc- One first 
constructs the real and symmetric matrix C defined by 



C^e-Y,lclI{B'e^Be), 



(9) 



which is diagonalizcd by the orthogonal matrix K such 
that D = KCK^, with Dxfj. — DxSxfi- The necessary 
transformation of the i?-matrix parameters is then given 

by i 



and 



E'x^Dx 



7c = ^7c 



(10) 



(11) 



It is straightforward to verify by substitution into 
Eqs. (PjItI) that these transformations leave U invariant. 



III. THE ALTERNATIVE 
PARAMETERIZATION 

A. Definition of tlie parameterization 

We begin by defining the real and symmetric matrix 



£ = e-J2lclIiSc~Be), 



and consider the eigenvalue equation 
Sat = Eitti 



(12) 



(13) 



where Ei is the eigenvalue and a.; is the corresponding 
eigenvector. Note that £ is implicitly dependent upon 
Ei through S'c, so the eigenvalue problem is nonlinear. 



We will assume for convenience that the eigenvectors are 
normalized so that aj ai — 1. 

Before proceeding further we would like to point 
out two important properties of this eigenvalue equa- 
tion: (1) The eigenvalues Ei are invariant if the Be are 
changed and the Ex and 7ac are changed according to 
Eqs. (UQjO). This result is easily shown by substituting 
Eqs. (@0) into Eqs. (|l|Q. (2) If = Sc{Ex), the ma- 
trix £ is diagonal for the energy Ex and hence Ex is an 
eigenvalue. For this choice of B^, the i2-matrix level en- 
ergy Ex is often taken to be the "observed resonance en- 
ergy" . This definition is particularly useful in the present 
context and we will thus adopt the Ei as the observed res- 
onance energies. The Ei also correspond exactly to the 
level energies found using boundary-condition constant 
transformations yielding Be = SAEx) such as described 
by Barker Q and Azuma et al. g. 

In addition one can define a new set of reduced width 
parameters 7^, 



via 



T 



(14) 



These new reduced width parameters are also invariant 
under changes in Be- When Be — Se{Ex) we have also 
7Ac = 7ac- The quantities Ei and 7,;c can be taken as an 
alternative parameterization of i2-matrix theory. We will 
derive below efficient methods to convert Ei and ^ie into 
the standard i?-matrix parameters Ex and 7ac, or to the 
collision matrix U. Also note that Ei and ^ie are equiv- 
alent to the "superscript (A)" parameters of Barker [||, 
and essentially equivalent to the "observed" i?-matrix 
parameters described by Angulo and Descouvemont |^ . 

Our Eq. ( p^ is closely related to the complex eigen- 
value equation introduced by Hale, Brown, and Jarmie 
to locate the poles of the collision matrix - in fact it is 
just the real part of their eigenvalue equation. For bound 
states our Ei are thus equivalent to the eigenvalues dis- 
cussed in Ref. [|j since Pe — 0. For these states we can 
also introduce the asymptotic normalization constant Cie 
which is given by H 



C. 



7' 

/ IC 



he \ dE I Ei 



(15) 



where ^c is the reduced mass. This quantity is simply re- 
lated to the pole residues described by Eq. (4) of Ref. [|| . 
For unbound states there appears to be no simple relation 
between Ex and 7ac and the pole parameters of Ref. [|j . 
One may however define the observed partial width of a 
level in terms of our parameters by 



r 



1 



l-ic He \ dE , 



(16) 



Ei 



see (LT, Eqs. XII.3.5 and XII.3.6). One should bear in 
mind however that there are many different definitions of 
observed resonance energies and widths in use; generally 
the differences between definitions are significant only for 
broad states. 
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B. Relation to standard parameters 

We will next show the method to convert E\ and 7ac 
to standard ii-matrix parameters. It is assumed that 
the eigenvalues are distinct, so that Ei ^ Ej provided 
i ^ j- Note that if this were not the case the levels with 
the same Ei could be combined into a single level. The 
eigenvectors of Eq. (|l^) are not orthogonal; using the 
eigenvalue equation with two different eigenvalues one 
finds 



aJ{Sj - £i)a., = (Ej - Ei)ajai 



(17) 



where £i is used to denote the matrix £ evaluated for the 
energy Ei. Using Eqs. (njjlj) with this result we obtain 



T 



E, 



(18) 



where Sic denotes the shift function Sc evaluated at Ei 
By similarly evaluating aj{£j + £i)ai, one finds that 



T 

ttj eai 



E, + E. 



'J T 

-a J ai 



(19) 



Sic + Sj 



-B, 



energy-independent matrices M and TV are completely 



determined by Ei, jic, and Be using Eqs. ( |2C| , |2l| ). The 
E\ can thus be determined by finding the eigenvalues of a 
generalized eigenvalue equation. If the matrix M is also 
positive definite then Eq. (|^) is known as the symmetric- 
definite eigenvalue problem and has A^a real eigenvalues 
(see Sec. 8.7 of Ref. M). The off-diagonal elements of 
M are ~ — J2c ^icljc^^ which is typically small com- 
pared to unity; M will be positive definite provided the 
7ic are not too large and the energy dependences of Sc 
are not too great. Further if M is not positive definite, 
the eigenvectors are not real and the transformation to 
standard i2-matrix parameters is not defined. We thus 
conclude that for physically reasonable 7ic, Ei, and Sic 
the matrix M will be positive definite; in practice we 
have found this condition to be easily fulfilled. Finally 
note that M is automatically positive definite for any 
given set of standard parameters since M = a^a. 

The eigenvectors of Eq. (24) bx can be arranged into 
a square matrix b which satisfies the relations 



and 



b^Mb = 1 



b^ Nb 



(25) 



(26) 



These results are summarized in the matrices M and N: 
{ 1 



T 

dj ai 



Ah, = 



Ei-Ei 



1= J 



(20) 



and 



T 

Uj eai 



(21) 



E^+J:clUS^c-Bc) 

- I El Sjc — Ej Si c 



B, 



^ = J 
i 7^ i 



Note that the construction of N requires the adoption of 
specific Be values. 

The eigenvectors of Eq. ( p^ ) can be arranged into a 
square matrix a such that Eq. (14) becomes 



7c 



T 

a 7c 



(22) 



The matrices M and N defined above can then be writ- 
ten as M = a^a and N — a^ea. From Eq. (^) the 
matrix e trivially satisfies the eigenvalue equation 



We therefore have b ^ a ^ and from Eq. (^2|) the stan- 
dard i2-matrix reduced widths are given by 



7c = b^l, 



(27) 



The simultaneous diagonalization of M and N thus pro- 
vides all of the standard ii-matrix parameters. Note that 
any Be can be chosen; the collision matrix U will be 
invariant provided the same Be are used in Eqs. ( |2l| ) 
and (^ or (^. The numerical solution of Eq. (^) is dis- 
cussed in Sec. 8.7.2 of Ref. [^; we have have utilized the 
LAPACK B routine dsygv. 



IV. FURTHER DEVELOPMENT 

It is fruitful to investigate alternative forms for the 
level matrix and the R matrix which allow the collision 
matrix to be calculated directly from the alternative pa- 
rameters. 



eux = E\u\. 



(23) 



A. The Alternative Level Matrix 



Upon substitution of u\ ~ ab\ and multiplying from the 
left by a^ this equation becomes 



Nbx = ExMbx 



(24) 



This eigenvalue equation holds the key for transform- 
ing from the Ei-"fic representation to the standard ii- 
matrix parameters Ex and jxc- The real, symmetric, and 



We define the alternative level matrix A implicitly via 
7JA7C = jJ,Aje. (28) 
In order for this relation to hold, wc must have 

aAa^ = A, (29) 
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or equivalently 



(30) 



where we have used Eq. (^). We can now substitute 
Eq. (0) for A^^ and again use Eq. (|2^ ) to obtain 



A ^ = a^ea — Ea^a 



A formula for the elements of Q "'^ in terms of the alter- 
native parameters can then be found using Eqs. ( ^3| | 



(39) 



(31) 



~ ~ Si^{E-Ej)-Sja{E-Ei 
^ic1]c Ei-Ei 



« 7^ j 



Using Eqs. (plJ3q) the collision matrix can now be written 



as 



N - EM -J2lclIiSc- B, + iPc). (32) U = 2ip"^0 \\ mPy^Rp^'^O^ +IO \ (40) 



The elements of this matrix can now be determined en- 
tirely from the alternative parameters with the aid of 

Eqs. mm- 




licljci.Sc + iPc) (33) 

C 



ic^ic 

~ ~ Sia(E-Ej)~Sj^{E~Ei) 
'^i<^'^]<: E^, 



Note that the boundary-condition constants Be have now 
canceled out - a not unexpected result since the alterna- 
tive parameters and the collision matrix are independent 
of Be- We can now express the collision matrix directly 
in terms of the alternative parameters using Eqs. (H,Pq) 



5c>c + 2i{Pc>Pcf'^llAic 



(34) 



B. The Alternative R Matrix 

The matrix R is an alternative to the standard R ma- 
trix and is defined implicitly via 



[1 - R{L - B)]-^R = (1 - iRPy^R, 



(35) 



where P is a purely diagonal matrix with elements Pc 
By comparison with Eqs. (^,^^ we must have 



(1 - \RP)-^R = 7^A7. 



(36) 



We proceed by assuming that R can be written in the 
form 



R = 7^Q7. 



(37) 



In the Appendix we describe a method to derive the level 
matrix form for the collision matrix [Eq. (H)] from the 
channel matrix form [Eq. (^]. This reasoning can also 
be applied to R and A. We find that in order to satisfy 
Eq. (^) we must have 



With the R matrix defined by Eqs. (|^,|9|) this equa- 
tion also gives U in terms of the alternative parameters 
without reference to the boundary condition constants. 



C. Relative merits of R and A 

The R matrix is more complicated than R and the 
calculation of U via Eq. ( pO| ) requires inverting a real 
-^A X N\ matrix in addition to a complex x matrix. 
When calculating U via the alternative level matrix one 
must invert a single complex N\ x N\ matrix - using the 
alternative R matrix approach may thus offer a modest 
computational advantage in comparison when N\ ^ N,.. 
Note however that if it is necessary to calculate U for 
several energies and N\ > Nc it will probably be more 
computationally efficient to diagonalize Eq. (|2J) once and 
then use the standard i2-matrix parameters in Eq. (^ to 
calculate U , as Eq. (||) only requires inverting a single 
complex Nc x Nc matrix. 

We would also like to point out that this alternative 
parameterization, using Ei and 7ic with Eq. (^) or (|4^), 
may be of formal interest since no arbitrary boundary 
condition constants are required, but the equations are 
mathematically equivalent to the standard ii-matrix ap- 
proach. The alternative parameters in fact correspond 
to eigenfunctions satisfying energy-dependent boundary 
conditions - the real part of the Kapur-Peierls or Siegert 
boundary conditions see (LT, Sec. IX. 2). 



V. SOLUTION OF THE NONLINEAR 
EIGENVALUE EQUATION 

The transformation from Eic and 7ic to the standard 
i2-matrix parameters E\ and 7ac can be carried out in 
a straightforward man ner using the methods discussed 
above in Subsec. IIIB. We will now discuss the inverse 



transformation, i.e. the solution of the nonlinear eigen- 
value problem Eq. (p^. At this point it is instructive 
to introduce a concrete example: in Table | we show a 
simple well-documented set of standard i?-matrix param- 
eters taken from Azuma et al. [^. 

We consider the linear eigenvalue equation 



Q ' = A-' +i')P^ 



(38) 



£{E)ai = EiCii 



(41) 



TABLE I: Standard Ji-matrix parameters from Table III of 
Ref. I which describe J" = 1 reactions in the ^^O sys- 
tem, and the alternative parameters derived from them. The 
channel labels a and 7 describe ^^C + q and ^^O + 7, re- 
spectively. The channel radius is 6.5 fm and the boundary 
condition constant is chosen so that the level shift vanishes 
for E = E\. The /3-decay feeding amplitudes B\ are equiva- 
lent to the quantities Aw'^^^Na^^'^ of Ref. [S]. 



A 


1 


2 


3 


£a (MeV) 


-0.0451 


2.845 


11.71 


7AC. (MeV^/") 


0.0793 


0.330 


1.017 


7A^ (MeV-^) 


8.76 X 10"^ 


-2.44 X 10"'' 


-2.82 X 10"® 


Ba 


1.194 


0.558 


-0.629 


£a (MeV) 


-0.0451 


2.400 


8.00 


7Aa (MeV'/'') 


0.0793 


0.471 


0.912 


7A^ (MeV-i) 


8.76 X 10"^ 


-3.20 X 10"" 


-2.50 X 10"® 


Ba 


1.194 


0.408 


-0.781 



where 5, the eigenvalues and eigenvectors hi de- 
pend upon on the energy parameter E. The solutions to 
the original nonlinear problem Eq. ( |l3|) thus occur when 
Ei{E) — E in which case E = Ei. From inspection of 
Eqs. (^,^2|) we see that the eigenvalues Ei also correspond 
to a set of standard i?-matrix level energies, transformed 
from the original parameter values to Be = Sc{E). 

We will next investigate how the Ei depend on E. 
Starting with 



E,, 



afSai 



differentiation with respect to E yields 



dE, 
dE 



—ea, + a,e— 



(42) 

(43) 
(44) 



+Ei i'^^cbi + of 

\dE ' dE 



da£\ 



Since by definition aj CLi 
and wc finally find 

dE, _ 
dE 



lwehave^d,-fdf ^=0 



■ 2 dSc 
' 'dE' 



(45) 



The energy derivative of the shift function is posi- 
tive for negative-energy channels, and is > for positive- 
energy channels for all cases we are aware of. This point 
is also discussed by (LT, p. 350); although a general proof 
of > is lacking it appears to always hold in practice 
and we will thus assume it is true here. Note that for any 
specific case it is a simple matter to verify this relation 
numerically. 

Since (7^0,^)^ is clearly > 0, we can utilize > to 



10 



> 







5 

E(MeV) 
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FIG. 1: The eigenvalue trajectories are shown by plotting 
as solid curves the eigenvalues Ei versus E; the dashed line 
corresponds to Ei — E . The eigenvalues of Eq. ( [l3| ) are given 
by the intersections between the solid curves and the dashed 
line. 



conclude from Eq. (45) that 



dE, 
dE 



< 0. 



(46) 



The eigenvalues Ei are thus monotonically non-increasing 
functions of E. The eigenvalue trajectories for the exam- 
ple parameters are shown in Fig. ^ where the expected 
behavior is seen. We also note the the eigenvalue trajec- 
tories avoid crossing one another for the reasons given by 
von Neumann and Wigncr [|o| . The avoided-crossing be- 
havior is most apparent when there are two nearby levels 
with very different reduced width amplitudes. 

The nonlinear eigenvalue problem Eq. (|l^) and the 
parametric eigenvalue problem Eq. ( ^ ) are also closely 
related to a well-studied question in linear algebra: the 
modification of a symmetric matrix with known eigen- 
values and eigenvectors by a positive-definite perturba- 
tion. This question is analyzed for the single-channel 
case in Sec. 8.5.3 of Golub and Van Loan [^| and for the 
multi-channel case by Arbenz, Gander, and Golub [ pT| . 
The perturbation bounds on the eigenvalues derived in 
ReL imply that the Ei remain finite provided the Sc 
are finite - thus the eigenvalue trajectories do not have 
real poles for \E\ < 00. 

From the lack of poles and the monotonic dependence 
< we can conclude that each eigenvalue trajectory 
intersects with the line Ei = E exactly once. These inter- 
sections are shown graphically for the example in Fig. |^. 
We thus have the important result that the non-linear 
eigenvalue problem Eq. (^3|) has a number of real eigen- 
values exactly equal to the number of i?-matrix levels. 
A similar type of nonlinear eigenvalue problem has been 
investigated by Rogers |jl^; it may be that the methods 
described in that paper could be used to develop further 
understanding of the present problem, e.g. to investi- 
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TABLE 11: Elements bij of the transformation matrix b cor- 
responding to the parameters of Table |. 



i \ 


1 


2 


3 


1 


1.000 


0.0373 


0.0446 


2 


0.000 


0.9781 


0.2281 


3 


0.000 


-0.1466 


0.9933 



gate inner products and/or the linear independence of 
the eigenvectors. 

The eigenvalues of Eq. ( |l3| ) satisfy the characteristic 
equation 

det(£: -E1)^0 (47) 

which can also be written as 

det [e-El- -yiS - B)j'^] = 0, (48) 

where is a purely diagonal matrix with elements Sc 
which depend upon E. Using the methods described in 
Rcf. one can show that 

det [e-El- j{S - B)j^] = det(e - El) (49) 
X det [l - 7^(e - El)-^j{S - B)] . 

The eigenvalues thus satisfy 

det(e - El) det [1 - R{S - B)] = 0, (50) 

which may be a computationally-efficient approach for 
determining the eigenvalues when N\ > Nc since the 
calculation of det(e — El) is trivial. Note that Eq. ( |50| ) 
is the multi-channel arbitrary-i3c generalization of the 
resonance condition given by Eq. (14) of Rcf. 0]. The 
eigenvalues also satisfy 



det [1 - R{S - B)] = 0, 



(51) 



but this equation has poles in addition to zeros, and if 
there is a level A with Be = Sc{E\) (at most one level 
can satisfy this condition) it does not produced a zero. 

Rather than finding the eigenvalues by directly solving 
the characteristic equation, we have applied the Rayleigh 
Quotient Iteration method described in Sec. 8.2.3 of 
Ref. j|] to Eq. (|l3|), as this procedure yields the eigen- 
vectors as well as eigenvalues. Starting values for the 
eigenvalues and eigenvectors were taken as Ei = Ei + 
J2c^ic[^ciEi) — Be] and aji — Sji. Due to the nonlinear 
nature of the problem, the matrix £ must be updated at 
each step of the iteration. These procedures were tested 
with several single- channel and multi-channel parameter 
sets, and were successful in finding all of the eigenvalues 
in every case. We cannot rule out however that some 
cases may require more carefully chosen starting values. 
Once the Ei and are found, the "fic can be calculated 
using Eq. ( [1^ ) . In Table || we show for the example case 
the alternative parameters determined from the standard 



i?-matrix parameters. Note that the alternative param- 
eters are exactly the same as the i?-matrix parameters 
given in the last column of Tabic III of Ref. |^ which 
have been transformed to satisf y Be — Sc{E\) for other 
levels. As discussed in Subsec. [II A this equality is re- 
quired due to our definition of the alternative parameters. 
In Table || we show the elements of the matrix b for the 
example parameters. Finally we would like to point out 
that the methods discussed in this section should be gen- 
erally useful for the extraction of resonance parameters 
from standard i?-matrix parameters. 



VI. APPLICATION TO 7 RAYS AND (3 DECAYS 

We will briefly discuss the application of the alternative 
parameterization to reactions involving 7 rays and (3 de- 
cays. Gamma-ray decay processes are generally treated 
with first-order perturbation theory in i2-matrix theory, 
which implies that 7-ray channels are excluded from the 
sum over channels when constructing A, A, M, or N. 
Assuming that external contributions can be ignored, the 
collision matrix elements connecting 7-ray channels (la- 
beled 7) and non-7-ray channels (labeled c) are given by 
(LT, Eq. XIII.3.9) 

- 2ir!,(P,P^)i/2^7Ac7M7^AM (52) 

= 2ir!,(p,p^) 1/27^^7^. 

In the long-wavelength approximation the penetration 
factor for 7 rays is given by P.^ — P^^^^ where £ is the 
multipolarity. The observed 7-ray widths are described 
by Eq. (|6|) , where 7-ray channels are excluded from the 
sum in the denominator. Using the same reasoning de- 



scribed in Subsec. IV A the alternative expression for the 
collision matrix elements can be obtained using the re- 
placement 



7c -477 = ijAj. 



7' 



(53) 



where the alternative 7-ray reduced width amplitudes are 
related to the standard parameters via 



77 = ^^77 



(54) 



If the external contributions to the matrix elements are 
included using the formalism of Barker and Kajino [ p^ , 
the expressions for the collision matrix elements and ob- 
served widths become more complicated. However these 
quantities can still be written in terms of the alternative 
parameters using the above equations, noting that the 
7A7 above are the internal 7-ray reduced width ampli- 
tudes. 

The extension of the alternative parameterization to 
the description of /3-delayed particle spectra is straight- 
forward. A multi-channel formula for the particle spec- 
trum is given by Eq. (7) of Barker and Warburton p4j; 



7 



note that additional parameters must now be intro- 
duced, the /3-decay feeding amphtudes B\x. It is con- 
venient to form column vectors Bx from the B\x , so that 
B\xAxfj,lfj,c can be writte n as j'^ABx- Again using 



the reasoning of Subscc. IV A we have 



7lABx 



(55) 



where the alternative feeding amplitudes Bx are related 
to the standard parameters via 



Bx=b^Bx. 



(56) 



Note also that if Be — Sc{Ex) we have Bxx = B\x- The 
/3-delayed particle spectrum can now be calculated di- 
rectly from the alternative parameters by using Eq. ( |55| ) 
in Eq. (7) of Ref. . One could also convert to standard 
i?-matrix parame ters u sing Eq. ( |56| ) and the methods dis- 
cussed in Subsec. IIIB, and then calculate the spectrum 
using standard i2-matrix formulas. 

In Tabic ^ we also show the standard and alternative 7- 
ray reduced amplitudes and /3-decay feeding amplitudes 
for the example case. 



VII. CONCLUSIONS 

We have presented an alternative formulation of R- 
matrix theory based on the parameters Ei and de- 
fined in Subsec. [II A . This parameterization is a gener- 



alization of the ideas presented by Angulo and Descouve- 
mont [0 . The new formulation is mathematically equiva- 
lent to the standard i?-matrix theory |3| but there are no 
boundary condition constants or level shifts. The new pa- 
rameters can be converted to standard i2-matrix param- 
eters by diagonalizing Eq. (p^, or be used to calculate 
the collision matrix directly using Eqs. (^4|) or (pO|). We 
have discussed the solution of the nonlinear eigenvalue 
problem Eq. (^) which is needed to convert standard 
i2-matrix parameters to the new parameterization. Fi- 
nally we have briefly discussed the application to 7 rays 
and f3 decays. 

We can envision at least two uses for this new formu- 
lation in the fitting of experimental data. One applica- 
tion is the generation of starting parameter values from 
an outside source of spectroscopic information such as a 
level compilation or shell-model calculation. These latter 
sources generally do not supply standard i?-matrix pa- 
rameters but rather resonance parameters without level 
shifts. In the past the methods to incorporate these types 
of information have not always been optimal (e.g. Be 
could be chosen to make the level shift vanish for a repre- 
sentative energy, but not for all energies simultaneously). 
Another application is to use the alternative parameters 
as the fit parameters. The calculations can be made di- 
rectly from the alternative parameters using the meth- 
ods discussed in Sec. IV, or by diagonalizing Eq. ( p^ ) to 
find the standard i2-matrix parameters. The latter op- 
tion may be preferable if Nx > Nc, if observables must 



be calculated for many different energies. It should be 
noted that in data-fitting applications the collision ma- 
trix must be calculated repeatedly for different energies, 
and the extra computational overhead required will be 
negligible in comparison. With the alternative parame- 
ters it is very easy to fix known information about level 
energies and partial widths for any number of levels. 
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APPENDIX A 

The equivalence of the two forms of the collision matrix 
given by Eqs. (|,|) is discussed in (LT, Sec. IX. 1). The 
derivation is reviewed here, utilizing the matrix notation 
introduced in Sec. II. The same procedure is useful for 
the derivation of the alternative R matrix as discussed 



in Subsec. IV B 



We define Lq = L B and note the quantity [1 — 
R{L — B)]^^ in Eq. (^ can be written as 

(l-RLo)-^^ (Al) 
[Lo - (io7^)(e - El)-\jLo)]-'Lo. 



A useful identity is given by 



{X + ZYZ'^)-^ =X 



(A2) 



-x-'ziY-' + z'^x-^zy^z'^x-^, 

which holds for any square and invertible matrices X and 
Y which need not be of the same dimension |^. With 
the aid of this identity we obtain 



(1 - RLo)-' = L^' - Lo\Loj^)x 



(A3) 



[~{e - El) + (7Lo)io"HW)] ' hLo)L^' Uo 



= l+7^(e-£;i-7io7^)"'7io (A4) 
= l+7^A7Lo, (A5) 

where in the last step we have used Eq. (|^) for the defi- 
nition of the level matrix A. 
We can then write 

(1 - RLo)-^R = (1 + 7^^7^^o) X (A6) 
7'^(e - Eiy^j 
= -y^ie- Eiy^-y (A7) 
+j^A[-A-^ + (e - El)]{e - Eiy^f 



8 



where we have substituted —A -'^ + (e — EV) for jLo'y'^ . which proves the equivalence of Eqs. (||,^). 
SimpUfying this expression we finally have 



(1 - RLo)^R = j'^Aj (A8) 
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